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Abstract 

In the zero mode approximation we solve exactly the equations of motion for lin- 
earized gravity in the Randall-Sundrum model with a non-standard distribution of 
I matter in the neighbourhood of the negative tension brane. It is shown that the form 

■ of this distribution can strongly affect the coupling of the radion to matter. We be- 

lieve that such a situation can arise in models with a realistic mechanisms of matter 



(N 



Q ' localization. 

o 

^ 1 Introduction 
m 

Q . The Kaluza-Klein hypothesis has been discussed in theoretical physics for three quarters of 

! a century. In accordance with this hypothesis, space-time may have extra dimensions, which 
are unobservable for certain reasons. The explanation of this unobservability, which was put 



O 



O ' forward in the original papers by Kaluza and Klein, implies that the extra dimensions are 



o 



X 



compactified and have a very small size of the order of the Planck's length Ipi = 1/Mpi. 

In 1983 Rubakov and Shaposhnikov put forward a new scenario for Kaluza-Klein theories, 
which was based on the idea of localization of fields on a domain wall PJ. They have also 
Qh' proposed an ansatz for multidimensional metric, which is compatible with this hypothesis 

In the last years there appeared indications that scenarios of this type can arise in the 
theory of strings [31 IH El E] (see Ref. [Zj for a review and references). In this approach 
our three spatial dimensions are supposed to be realized as a three-dimensional hypersurface 
I embedded into a mutidimensional space-time. Such hypersurfaces are called 3-branes, or 
just branes. The main goal of such scenarios was to find a solution to the hierarchy problem. 
It was solved either due to the sufficiently large characteristic size of extra dimensions jH], 
or due to exponential warp factor appearing in the metric [Oj. In both approaches gravity in 
multidimensional space-time becomes "strong" not at the energies of the order of 10^^ GeV, 
but at much lower energies, maybe of the order of 1 10 TeV. An attractive feature of 
these models is that they predict new effects which can be observed at the coming collider 
experiments. 

In paper [H] an exact solution for a system of two branes interacting with gravity in a 
five-dimensional space-time E was found. This model is called the Randall-Sundrum model 
(usually abbreviated as RSI model), and it is widely discussed in the literature (see Refs. 
[TUl ITT] for reviews and references). Let us denote the coordinates by {x*^} = {x^,y}, 
M = 0, 1, 2, 3, 4, /i = 0, 1, 2, 3, the coordinate x'^ = y parameterizing the fifth dimension. It 
forms the orbifold /Z^, which is realised as the circle of the circumference 2R with points 
y and —y identified. Correspondingly, we have the usual periodicity condition in space-time 
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E, which identifies points {x,y) and {x,y + 2nR), and the metric qmn satisfies the orbifold 
symmetry conditions 

9,,u{x,~y) = gf,u{x,y), (1) 
gf^iix, -y) = -gf^iix,y), 
5-44(3;, -y) = guix,y). 

The branes are located at the fixed points of the orbifold, y = and y = R. 
The action of the model is 

S = Sg + Si + S2, (2) 

where Sg, Si and 5*2 are given by 



^ I671G 



^ / (R-A)^ d'xdy, (3) 

51 = Vi ^/^5{y)d'^xdy, 

Je 

52 = V2 [ ^5{y - R)d^xdy. 

Je 



Here g^j_y is the induced metric on the branes and the subscripts 1 and 2 label the branes. 
We also note that the signature of the metric guN is chosen to be (—,+,+,+,+). 
The Randall- Sundrum solution for the metric is given by 

ds^ = gMNdx'^ dx^ = '-yf^ydx^dx" + rfy^, (4) 

where 7^1, = e^'^'^^-'r^^i,, 77^,^ is the Minkowski metric and the function cr{y) = —k\y\ in the 
interval —R < y < R. The parameter k is positive and has the dimension of mass, the 
parameters A and Vi,2 are related to it as follows: 

A = -I2P, Vi = -V2 - 
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We see that brane 1 has a positive energy density, whereas brane 2 has a negative one. The 
function a has the properties 

d,a = -k stgn{y), ^ = -2k{5iy) - 5{y - R)) = -2kl (5) 

Here and in the sequel 84 = -^. 

We denote k = V IQnG, where G is the five- dimensional gravitational constant, and 
parameterize the metric gMN as 

guN = 7mn + f^hMNi (6) 

Hmn being the metric fiuctuations. Substituting this parameterization into © and retaining 
the terms of the zero order in k, we get the second variation action of this model ^21- It is 
invariant under the gauge transformations 

h'uNi.^^ y) = hjvmix, y) - iyuiN^x, y) + WNiu^x, y)), (7) 
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where Vm is the covariant derivative with respect to the background metric 7m7v, and the 
functions C,N{x,y) satisfy the orbifold symmetry conditions 

e{x,-y) = e{x,y), (8) 
^^{x,-y) = -^^{x,y). 

With the help of these gauge transformations we can impose the gauge 

^4 = 0) ^44 = hii{x) = 0(a;), (9) 

which will be called the unitary gauge (see f21)- We would like to emphasize once again 
that the branes remain straight in this gauge, i.e. we do not use the bent-brane formulation, 
which allegedly destroys the structure of the model, i.e. the orbifold symmetry under the 
reflection —y ^ y (this problem was discussed in jl3j). 

In the linear approximation interaction with matter looks like 

^ I h''''{x,y)TMNV^d'xdy, (10) 

where Tmn is the energy-momentum tensor of the matter: 

6L 

Thus, we are considering perturbations about the background, for which Tmn = 0. 

We will examine the /i^i^-components of the metric fluctuations huN, since all four- 
dimensional physical effects can be described in terms of this field. Obviously, the unitary 
gauge conditions © do not fix the gauge of this field. In fact, after imposing these gauge 
conditions there remain gauge transformations of the form 

= e^"^,{x), (11) 

which change the longitudinal components of the field h^y. Nevertheless, it turns out that 
it is convenient to solve the equations of motion for linearized gravity in the unitary gauge 
and then to choose an appropriate gauge in our four- dimensional world on the brane. We 
will use the de Bonder gauge for the field h^i, on the brane, which corresponds to the choice 
of harmonic coordinates. 



2 Equations of motion 

The equations of motion for different components of the metric fluctuations in the unitary 
gauge take the form (see [T^ ITH]): 
1) /ii/-component 

\ {dpdPh^y - d.d^hpy - dydPhp^ + - 2A;2 V + U^djh + (12) 

^ v 
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2) /i4-component, 

diidf^h - d^h^^) - 3d4adf,(j) = -kT^4, (13) 

which plays the role of a constraint, 

3) 44-component 



^(a^a^ V - d^d^^h) - h^adji + 6A;V = -'^T,,, (14) 



with Tmn being the energy- momentum tensor of the matter and h = In what follows, 

we will also use an auxiliary equation, which is obtained by multiplying the equation for 44- 
component by 2 and subtracting it from the contracted equation for /iz/-component. This 
equation contains h and only and has the form: 

2 + 2diad4h - Sk'^cj) + 8k(j)5 + d^d^cp = ^ (T„^ - 2Tu) ■ (15) 

For example, if Tmn = 0, the physical degrees of freedom of the model can be extracted 
by the substitution [12] 

hf,u = bf,u + If^uio- ~ + ^ (^^ " c + ^ + 1^"^'") d^ducj). (16) 

with c = ^2kR_i ■ It turns out that the field b^^{x^,y) describes the massless graviton and 
massive Kaluza-Klein spin-2 fields, whereas (f){x) describes a scalar field called the radion. 

However, the situation is rather different, when there is matter on the branes. The cases, 
in which matter is located on the branes, i.e. the energy- momentum tensor is of the form 
Tf_iu = tf^^{x)S{y) or T^,^ = tf^^{x)6{y — R) (and T^4 = 0, T44 = 0) were discussed in detail 
in ^j. It is well known that if we live on the negative tension brane (at y = R), 
the contribution of the radion is e^^^ times stronger than the contribution of the massless 
graviton ^1 . It means that in the case of the massless radion scalar gravity is realized on 
brane 2, and it is necessary to have a mechanism for generating the radion mass, for example, 
the Goldberger-Wise mechanism fT^ to make gravity in the zero mode approximation tensor, 
i.e. to suppress the contribution of the scalar radion component. But we will show below, 
that this is not the only solution for this problem. 



3 A simple example 

In the original formulation of the Randall-Sundrum model the mechanism of localization of 
fields is not taken into account. The branes are treated as infinitely thin objects, and it is 
assumed that energy- momentum tensor has 6{y)-like profile in the extra coordinate (matter is 
located on the brane only). Nevertheless in some papers (see, for example, jT7|ITHllT^I^ I^) 
considering models with extra dimensions, energy-momentum tensors with dependence on 
the extra coordinates were used for constructing non-trivial background metric. One can 
also recall models with universal extra dimensions, in which physical fields can propagate 
in the extra dimension [22j, and fat brane scenarios [211 1211 • Thus, the existence of energy- 
momentum tensors of this type seems to be quite reasonable. Here we consider some simple 
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examples, which are sufficient to show attractive features of the non-delta-hke locahzation 
of matter. 

Let us choose as an example the following form of the energy-momentum tensor 

= t^^{x)^-, (17) 
T^4 = 0, (18) 
T44 = t,,{x)Y''e-'''^{a + e''^c-ec), (19) 

where e is some arbitrary constant, which will be defined later. Here t^j_y{x) is the energy- 
momentum tensor of matter on the brane, it depends on the four-dimensional coordinates x 
only. One can see that the matter is localized near the negative tension brane at y = R. We 
would like to note that the energy- momentum tensor of the form (jl7p . (jl8p and (jl9|) satisfies 
the covariant energy conservation law, which has the form 

V^Tmtv = 0, (20) 

and the function of localization in (jl7|) is normalized to unity. Constructions similar to 
that used in this paper, but in the case of six dimensions, were utilized in [23 12^1 • The 
explanation of such constructions depends on the method of localization, and this issue will 
not be discussed in this paper. 

The substitution, which allows one to decouple the equations (fT^ . (fTSj) . (fT^ . (fT3j) with 
Tmn of the form (II7|), (HHD, dH) looks like 

where h^j_i, was defined in and u^j_p{x^,y) describes the massless graviton and massive 
Kaluza-Klein spin-2 fields. Substituting (fTT)|) and (PT|) into equations (fT!?jl . (fT^ and (fTHjl . we 
get: 

94(6-2^^(9%, -(9^m)) = 0, (22) 

e-^^id^'d'H^, - Uu) - 394^x94 (e-2"M) + e-^-^—n^ - (23) 
e ^ ' R ^ ' 

^4 (6^-94 (e-2-M)) + in^e-^'^ (1 - 2 (a + ce''^^)) = (24) 

= |^te-2-(l-2(a + ce2'=^)), 

where m = rf^Up^y, t = r]'^'^t^,y, d'^ = rj^'^djy, □ = 7]^''d^dy and 

c = c(l-ee-2'=^) 
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Let us consider the Fourier expansion of all terms of equation (j24j) with respect to co- 
ordinate y. Since the term with the derivative has no zero mode, this equation implies 
that 

□0 = (25) 

54(e^2"M) = 0. (26) 

Now let us consider /xi/-equation. It is well known that the field Uy,y in the presence of matter 
is a combination of zero and massive modes, whose eigenfunctions are orthogonal |12|. In 
particular, the zero mode can be represented as = e^'^ a^y, where a^^y depends on x only. 
It also means that with the help of the residual gauge transformations (fTT|) it is possible to 
impose the gauge condition 

d^u^y = 0, (27) 
«7 = 0, 

on the massive modes u'^^, and the de Bonder gauge condition on the zero mode 

d" (a^u - ^Vf^uO^ = 0. (28) 

Having imposed this gauge, we are still left with residual gauge transformation 

= e'^'e^ix), ne^ = 0. (29) 

The gauge transformations with satisfying these conditions are important for determining 
the number of degrees of freedom of the massless graviton. It follows from equation 
that 

^ 2A;R 2 

Da = 2 — t. (30) 

R 

Substituting (fTHjl and (pT|l into equation (fT^ with condition and passing to gauge (P7j) . 
we get 

lo{a,. - + ^e-2-nn;r, + \d^d^u^^, - 2^^^™ - d^d^au^^, = (31) 



2R 

We are going to calculate the equations of motion in the zero mode approximation. Thus, 
we have to find an equation for the field If we multiply equation (jHT|) by e^°", integrate 
it over y and take into account the orthonormality condition for the wave functions of the 
modes [T2j, we get 

1 kc^e^^^ 
°(aMi^ - i^Vf^va) = -2 — - — tf,^. (32) 

The equation for massive modes takes the form 

ie-'^"D<„ + is.aiu™ - ik'w;:, - ftSiKl = (33) 



Now we are ready to find equations of motion for the zero mode part of /i^^,. Since the 
matter is locahzed near the brane, we will calculate the equation for /i^;y|y=_R. It is easy to 
see that h^y\y=ji (fTBj) . (PT|) does not satisfy the de Bonder gauge condition. The residual 
gauge transformations (j^^ are not sufficient to pass to this gauge. But since we consider 
only the effective theory on brane 2 (at y = i?), we can drop the condition De^ = 0, which 
fixes the gauge for the field h^y. Then we can pass to the de Bonder gauge condition for the 
field hp^p on brane 2 with the help of these gauge functions. Making these transformations 
(analogously to what was made in jlllllS]), we get in the zero mode approximation 



(^V - \y=n = (34) 



An important point is that these equations are written in the coordinates {x^}, which 
are Galilean on brane 1 (not on brane 2) and are inappropriate for studying physical effects 
on brane 2 (we recall that coordinates are called Galilean, if g^,^ = 

diag{-l,lAA) EH)- It 

is necessary to pass to Galilean coordinates on brane 2 to get a correct result. This problem 
was discussed in detail in papers ITil IT^ . In Galilean coordinates on brane 2 equation 
dMl) looks like 



(^V - ^Vf^uhj \y=R = (35) 



e-'"'' ( - lv,ua] - 2{-kR - c) ^^^^ 



and equations (j2SI), take the form 



□0 = (36) 
^{(^t^y - i^Vi.ua) = -2 — - — tf,^. (37) 



Thus, we can get 



° I - 2^t^'^h \ \y=R = (V - 3 yii^u - jtj, (38) 

where t^u is the energy-momentum tensor of matter in the coordinates, which are Galilean 
ai y = R. One can see, that this equation coincides with equation for the fluctuations of 
metric in the linearized Brans-Bicke theory, which looks like 

□ (Sg,. - lv,uSg^ = -16nG (^V " {v,. ' ^) , (39) 

where u) is the BB-parameter and G is the gravitational constant. Thus, comparing equations 
dSHl) and dSni), we get 

n„2 2kR 

G = G-^, (40) 
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.= 5(1^. (41) 

It is easy to see, that in the the Randall- Sundrum model with the (5-function-like localization 
of matter on the brane there is a factor e^^^ instead of e in ()38|) (compare with the analogous 
formula in jl4t [T3] ) . 

From the recent experimental data (see, for example, |2H|) we know, that uj > 3500. It 
means, that 

e<4,3-10~^ (42) 
With this value of e there are already no problems with the radion field. 

4 Localization of matter 

In (fT7j) the factor e~^'^, describing distribution of matter in the extra dimension, was utilized. 
However one can say that matter is spread in the whole bulk rather than localized on the 
brane (in analogy to the massless graviton whose wave function is ~ e^'^), and it is not 
evident, why we make all calculations for h^^, aX y = R. It seems that it is more reasonable 
to consider matter, which is confined to brane 2 much stronger, since the description of 
gravity on brane 2 by h^y\y=R is more justified in this case. Let us consider a "toy model" 
with energy-momentum tensor of the form 

T =t (x) ^ e~^"(^+^^ f43) 

J-fiu ^l^i^y-^) ^2kR+2NkR _ ' 

where the function of localization is normalized to unity as well. Other components of the 
five-dimensional energy-momentum tensor are chosen to be 

T,, = 0, (44) 

T44 = t{x) {k + kN) e-'^ (^^- (45) 
cee^NkR ^ 



Q2kR+2NkR _ ^ 2Nk (^e^kR+2NkR _ -[^ 



^-2Nc 



Such a complicated form of (j45|) is caused by the requirements to satisfy the energy conser- 
vation law ()20|) and to get a weak coupling constant for the radion. 

The substitution, which allows one to decouple and solve equations (fT^. and nv- 
equation with PHj) . (jl^ and (jlSj) . has the following form: 

2NkR 

2cj~2Na 




ANk {Nk + 2k) e 
Analogously to what was made in Section 3, we can get 



kt {k + kN) e^^^^ 

~ 3 (^^2kR+2NkR _ ]^-) ^' (47) 
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1 kk{k + Nk) {1 - e-^^''^) 

Using ()46|) and passing to Galilean coordinates at y = R, in the zero mode approximation 
we get 



° (^V - lv^.uh^ \y=R = (49) 



kc{l + N){l-e-^^'''^) f 2NkRt ( d^d, , 



NR (1 - e-2fc«-2^^«) V 3 (1 - e-2^^-f^) V □ 



One can see that in the limit — > equation ()49|) passes into ()38p . The four- dimensional 
gravitational constant has the form 

AT/? (1 - e-2fc«-2^'=«) ^ ^ 

(compare with (jHH))- For a large we get 

G^G^. (51) 

A choice of relatively small e makes the contribution of the radion to be not in contradiction 
with the present-day experimental data. 



5 Conclusion and final remarks 

There may arise a question about generation of the new hierarchy instead of the one solved 
in the Randall- Sundrum model. Actually, the factor e in (jl9|) is much smaller than the factor 
g2fc_R jj-^ (ji9j). But for the energy-momentum tensor on brane 2 (at ?/ = R) we get 

g2fc_Rg 

T \ r, — f ~ hf 

ri 

Tii\y=R = t—— (l - 4.3 ■ 10'^) ^ t k'^R (l - 4.3 ■ 10"^) , 
R 

and the correction ~ e is not so small in comparison with the main value ~ 1 (analogous 
calculations can be made in the case (jIHj) and (jlSj) too). It seems to be quite reasonable. 
One should take into account that the energy- momentum tensors of the form ()17|) . (|T8|) . 
(fTIHl and pHjl . (jllj), PKjl do not correspond to any real field action, and must be interpreted 
as "toy models". There also arises the question about the origin of such correction ~ e. 
Explicitly the term ~ e~^'^ in (fT^ and corresponds to the homogeneous solution of 
equation (j^Uj) (i.e. to the case of T^^, = 0). Nevertheless such term can be taken into 
consideration, since its existence does not evidently contradict our assumptions. Surely it 
would be better to derive similar form of the energy-momentum tensor in a more natural 
way. For example, simple examples with scalar and gauge fields in the bulk in the Randall- 
Sundrum model were discussed in papers I^HEHI- Wave functions of these fields have their 
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maximum values on brane 2, and these solutions can be treated as a basis for our choice of 
the energy-momentum tensor. But even in the simplest case of the scalar field only, there 
arises a nonzero T^4-component of the energy-momentum tensor, which does not satisfy the 
condition ()18|) . Moreover, the wave functions are different not only for different fields, but 
even for different modes. Thus, formulas (fTTj) . (fTHj) . (fT^ and (jlHj) . (jH)), must be used 
as "toy models" for classical objects, but not for the fields (for example, one can choose 
too ~ m6{x), tio = 0, Uj = for a static point-like source). 

In this paper we have discussed only one aspect of the radion field problem. Indeed, 
though the coupling constant of the radion to matter can be small enough to fit the experi- 
mental data, the size of extra dimension in the RSI model should be stabilized in any case. 
It is well known that in this model gravity alone cannot stabilize the size of extra dimension 
due to the Casimir effect But if matter is not localized exactly on the brane, i.e. it has 
wave functions in the extra dimension, it is possible that the contribution of matter fields 
to the Casimir force between two branes could stabilize the size of the extra dimension. 
Moreover, even in the case of classical stabilization with the help of scalar field (such as in 
|32j ) the coupling constant of the massive radion to matter can vary essentially depending 
on the profile of matter in the extra dimension. 

There is another question, which may arise in connection with our choice for matter 
distribution: why do we choose non-delta-like profile for matter, whereas the branes have 
delta-like profile? Of course, such profile of the branes is an idealization. But one can recall 
kink-like solution found in [2H], which describes a delta-like brane in the thin brane limit. 
Thus in the case under consideration the RSI setup can be regarded as a phenomenological 
model for describing the branes in the case when the thickness of branes is much smaller 
than the effective width of the wave functions of matter. Another advantage of the RSI 
background is its relative simplicity, which allowed us to decouple and solve exactly equations 
of motion for the radion field and tensor modes at least in the zero mode approximation and 
to show possible effects which can be produced by a non-standard matter distribution. In a 
more realistic model the non-zero width of the branes should be definitely taken into account, 
but it is evident that effect analogous to that described above should exist in this case too. 

It should be also noted that the longitudinal parts of substitutions (fT^ . (PT|) and (jlUj) 
contain some terms, which can lead (but it is not necessarily so) to the strong coupling 
effect, analogously to what happens in the DGP model jSHES]- These terms are pure gauge 
from the four-dimensional point of view from the brane and are not dangerous in the linear 
order. Nevertheless it is necessary to consider the non-linear corrections to the equations of 
motion to discover the contribution of these longitudinal terms. But this problem deserves 
an additional detailed investigation. 
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